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Abstract 

Evolutionary structure searches are used to predict a new class of compounds in the lithium-rich 
region of the lithium/hydrogen phase diagram under pressure. First principles computations show 
that Li^H, 4 < m < 9, are stabihzed with respect to LiH and Li between 50-100 GPa. The building 
block of all of the lithium subhydrides is an LigH cluster, which can be thought of as a superalkali. 
The geometries and electronic structures of these phases is analogous to that of the well-known 
alkali metal suboxides. 
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Considerable attention has been given towards the pressure-induced synthesis of 
hydrogen-rich systems with unusual stoichiometrics due to their potential import in en- 
ergy applications, and because they may provide routes towards hydrogen's metallization 
and concomitant superconductivity. [Ij A number of intriguing systems which exhibit strong 
intermolecular interactions under pressure, including the van der Waals solids Xe-H2,[2] 
CH4-H2,[3j SiH4-H2,[4i |5j and H2S-H2,^ have been prepared. Pressure-induced hydrogen 
uptake has been observed in metals with low hydrogen solubility yielding RhH2,[7j Re2H,[8] 
PtH.[9l-[TT] or noble metal hydrides. [12j Yet, not much work has been devoted to studies of 
compounds which are hydrogen poor. 

Recently, we have shown that under pressure the alkali metal poly hydrides, MH^ with 
n > 1, become stable with respect to decomposition into MH and H2 at 100, 25, and 2 GPa 
for lithium, [13j sodium, [T4], and rubidium, [15j respectively. LiH^ is particularly interesting 
because the most stable structure above 150 GPa, LiHe, was found to be metallic already at 
1 atm.[13] Elemental lithium also exhibits unusual behavior: when compressed its melting 
temperature initially decreases so that between 40-60 GPa it has the lowest melting point 
among the elemental metals. [16] These reasons prompted us to computationally explore 
the alkali-metal/ hydrogen phase diagram further, specifically its lithium-rich region. Even 
though the structures we find, some of which are illustrated in Figure [T| show no evidence 
for significant metallicity or low-temperature melting behavior, they hint that pressure may 
be used to synthesize a new class of compounds whose atomic and electronic structures are 
remarkably similar to the well-known rubidium and cesium suboxides. [171 [IE] Because of 
the similarities described below, we refer to these phases as lithium subhydrides. 

The open-source evolutionary algorithm (EA) XtalOpt [20j has been employed to find 
low enthalpy structures of Li^H with 2<m<9uptoP = 100 GPa.[21J The calculated 
enthalpies of formation, AHp^ of the most favorable structures for a given stoichiometry 
are provided in Figure [2](a). The AH^ values are computed from the enthalpy diflFerences 
between the right- and left- hand sides of the reaction shown at the top of Figure [2} quan- 
tum and temperature eflFects are not included in this plot. In this pressure range (between 
50-100 GPa) a number of Li^H extended systems with m = 4 — 9 are found to become (en- 
thalpically) stable with respect to decomposition into Li and LiH, their would-be precursors 
in experiment. Moreover, the formation of LiH from lithium and molecular hydrogen is pre- 
ferred up to at least 350 GPa [13j. Our computations predict that two different LisH phases. 
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FIG. 1: (a) Extended structures [19] of the two isoenthalpic LisH structures found between 50- 
100 GPa and (b) the chain-hke building blocks common to both systems. These chains are made 
up of (c) LigH clusters which are fused at the faces as illustrated in (d) by the Li^ atoms. 



one with Cmc2i (space group #36) and one with Abm2 (space group #39) symmetry (see 
Figure [l]( a)), are particularly stable in the sense that they have the most negative enthalpies 
of formation between 50-100 GPa, a pressure range which coincides with the metal-to- 
semiconductor transition [22j and the low melting point [16j observed in pure lithium. Both 
structures have enthalpies within 3 meV/atom of each other within the entire pressure range 
studied. Phonon calculations confirm that both are mechanically stable within the harmonic 
approximation at 90 GPa and have the same zero-point energy. Beyond 180 GPa a Li2H 
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FIG. 2: AHp for the reaction LiH+ ((m — l)Li) Li^H versus pressure for various integer values 
of m. The enthalpy for Li is calculated for the most stable structures from Ref. [23j and the 
enthalpy vs. pressure plot of the Li phases we considered is shown below. In this pressure range 
LiH is found in the NaCl structure. 

stoichiometry which is best described as a CsCl structure with missing hydrogen atoms has 
a lower AHjp^ but since this is also a region where lithium polyhydrides are stabilized. [T3] 
we will focus our discussion on m > 3. 

Structurally, all of the Li^H with 3 < m < 9 are related to each other through common 
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motifs: LisH cluster-like building blocks. The hydrogen atom lies near the center of these 
clusters and the lithium atoms make up the vertices of their enclosing polyhedra, the faces 
of one such polyhedron are colored red in Figure |T[c). The most common cluster's geometric 
structure, which resembles a previously studied LigH cluster, [24j is shown in Figure [l](c), as 
extracted directly from the LisH extended system. It is present in all of the structures except 
for m = 4.[25j This C^^symmetry polyhedron has 11 faces (10 triangular and 1 rectangular) 
and is best described as a distorted bicapped trigonal antiprism. Note that two opposing 
Li-Li distances on its "rectangular" base are distorted such that there is no C2 axis in this 
cluster. 

The LigH unit is also qualitatively similar to the most stable LigB isomer found in an 
extensive investigation of Li^B clusters. [26j A [LieBj^ cluster was found to be particularly 
stable relative to the rest and was noted to correspond to a closed electronic shell (magic 
number) within the phenomenological shell model, |27j as does [LigHj^. We found this impor- 
tant because the building blocks of the suboxides, Rb902 and CsnOa, have been compared 
to giant alkali metal (superalkali) atoms which may combine with excess metal atoms yield- 
ing a type of intermetallic compound. fT8] The neutral building block of the subhydrides, 
LigH, has the correct electron count (8 + 1) so that it can be thought of as a superalkali 
atom, a cluster with one more electron than a magic number (eight in this case). Moreover, 
interstitial lithium and/or hydrogen atoms (atoms which do not belong to a LigH building 
block) are observed in Li^H with m = (3, 6 — 9) (see Figs.|3](b) and (c)). Thus, the extended 
systems predicted in this work have both structural and electronic similarities to the subox- 
ides. Interestingly rare-earth halide cluster complexes almost always possess an endohedral 
atom as well. [2^8] 

In LisH, the LigH clusters are fused at the faces in a chain-like arrangement, as illustrated 
in Figure [3](a) ; the two LisH phases only differ in the relative stacking of the chains (see 
Figure [l](a,b)). LiaH and LieH are constructed in a similar fashion, except in LieH the 
polyhedra are fused via edges and vertices (see Figure [3](b)). Similarly, Li4H can also be 
thought of as a network of edge- and vertex-sharing clusters. Li^H with m=7-9 are made 
up of pairs of face-sharing polyhedra, as illustrated for LiyH in Figure [3](c). The range of 
Li^H structures we considered (all of which are predicted to be thermodynamically unstable 
with respect to decomposition into LisH and Li/LiH) can all be represented as arrangemnts 
of LigH building blocks with lithium atoms possibly in the interstital regions. 
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FIG. 3: The underlying chain-like polyhedral building blocks of (a) Cmc2i-Li^}i and (b) LieH 
at 90 GPa. The fused polyhedron pairs in LiyH at 90 GPa are shown in (c). 

Here lies the most striking similarity to the alkali metal suboxides: a structural one. 
The suboxides' basic building blocks, Rb902 and CS11O3, are formed from two RbeO or 
three CseO octahedra; they are similar to the fused polyhedron pairs shown in Figure [3](c) 
and are the only constituents of the simplest suboxides. Different stoichiometrics, such as 
Rb60=Rb902-Rb3 and Cs70=Csii03-Csio, are made up of the same basic units with excess 
Rb or Cs atoms interspersed between them, like LiyH discussed above. The underlying 
octahedra can also be fused such that they form chains in more complex subnitrides or 
suboxoindates, such as NaBaaN [17j or Cs9ln09,|29j just like LisH. 

In analogy to the suboxides and subnitrides, one can assume that because of the large 
electronegativity difference between lithium and hydrogen, the lithium subhydrides can be 
thought of as (Li+)^(H~)-(e~)^_i, exhibiting ionic bonding on the inside and metallic bond- 



6 



ing on the outside of the fused LisH clusters yielding an overall metallic character. The cal- 
culated Bader charges of -0.96 on the hydrides and +0.68 on the lithium atoms agree with 
the ionic bonding, and the ubiquitous presence of the LisH building blocks suggests strong 
intra-cluster atomic interactions. The DOS of the two LisH phases at 60 and 90 GPa show 
these are semimetals or are weakly metallic, see Figure [4](a). Furthermore, the existence of 
two iso^enthalpic LisH phases suggests that, like the suboxides, the bonding between the 
chains is weak. 




Energy / eV 

FIG. 4: (a) The total electronic densities of states (DOS) for Cmc2i LisH at 90 GPa (solid line) 
and 60 GPa (dashed line). The DOS plots for Abm2 LisH are qualitatively similar and are given 
in the SI. The left and right plots have different scales, (b) The site-projected densities of states 
(PDOS) onto the Lia and Li^ atoms in Cmc2i-Li^}i at 90 GPa. The Fermi energy, E^, is set to 
zero. 

Under pressure the valence bands of LisH narrow slightly, the semi-core (sub-valence) 
bands do not change much, and the core Li Is bands broaden, much like in elemental 
lithium. Also like lithium, the density of states near the Fermi energy, E^^ has mostly Li p 
character. However, at 60 GPa LisH has a lower DOS at Ep than Li c/16 and at 90 GPa 
the subhydrides remain semimetallic but Li-o(740 is a semi-conductor. The band^structure 
of Cmc2i Li^}l in Figure IsFa) shows that along the high-symmetry lines two conical-like 
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bands, with a linear dispersion, cross E^. In Abm2~Li^]i the band crossing occurs shghtly 
above Ep^ so this phase is somewhat more metaUic. 

The bands lying between -9 and -6 eV are primarily hydridic, but some lithium character 
is observed. The Li^, atoms, which belong to a face that is shared between two adjacent LisH 
clusters, contribute more than the Li^ to the low energy portion of the sub- valence bands, as 
shown in the PDOS in Figure [4](b) . In order to investigate this further, a Wannier function 
(WF) was constructed which describes these bands almost perfectly (see the SI) using the 
A^MTO method. [30] This WF, shown in Figure [5](b) , consists of a sphere centered on the 
central H~ with pAike tails on the surrounding lithium atoms, indicative of hybridization 
between the two. Interestingly, this WF is reminiscent of a Is superatom orbital in a [LigH]^ 
cluster with the same geometry as the building block used to construct the LisH phase (see 
the SI). It also shows that the bonding within the cluster is primarily ionic in nature. 




FIG. 5: (a) Band structure of Cmc2i-Li^}i at 90 GPa, (b) the Wannier Function which yields the 
two sets of bands between -9 and -6 eV (isovalue= zb 0.04 a^^^^, where is the Bohr radius), and 
(c) an ELF=0.87 isosurface. In (c), the two characteristic maxima in the interstitial region of the 
lattice are labelled Ml and M2 and are shown in greater detail in the SI. 

We note that even though the WF which describes the sub-valence band of LisH in Figure 
[5](b) is reminiscent of a superatom 5-orbital, the valence DOS cannot easily be understood 
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using such a picture. To extend the analogy, the remaining occupied bands would need 
to be described by superatom p-like states. Chain-like cluster arrangments of LigH, such 
as LiisHa, do qualitatively behave as such (see the SI), but we will show below that the 
comparison of the Li^H phases with the suboxides (and high-pressure lithium) is perhaps 
a stronger one. 

As expected, the electron localization function (ELF) is high around the hydridic atoms. 
There are, however, two characteristic pockets of increased electron density, labelled Ml 
and M2, with peak ELF values near 0.91 between the chains (see Figure [5](c)) that can be 
attributed to the metallic bands. The ELF profile of the Abm2-Li^}i structure is similar 
and is consistent with the localization of valence electron density in the interstitial regions 
noted in pure compressed lithium. [311 [32 J The integrated charge densities within the M1/M2 
basins (the charge density maxima identified from the Bader analysis) are 0.74/0.84 electrons 
in both LisH structures, notably smaller than the 1-2 electron occupancies in the recently 
proposed C2ch-AQ structure for the pure hi-oCAQ phase found experimentally at 90 GPa; 
this is suggestive that LisH is more metallic. [23] as expected. However, the ELF between 
two M2 basins reaches values as high as 0.74, and the ELF between Ml and M2 does not 
exceed 0.68. The diflFerence between these ELF values and the basins' peak ELF values of 
0.91 may be suggestive of poor conductivity. [231 [33 J also as expected. 

The localization of metallic electrons in the areas between the chains are also character- 
istic of the suboxides. [18j It is interesting to note that the Ml site shown in Figure [5](c) is 
actually in the same local chemical environment as the hydrogen atoms. That is, the Li 
atoms around the site form the same exohedral structure as they do in the LigH cluster 
discussed above. Given that the 0.74 electrons assigned to the Ml site is comparable to the 
-0.96 charge assigned to the hydrides, this provides another trait (beyond the 1-dimensional 
behavior of the hydridic bands) which is unique to LisH among the Li^H we considered: it 
allows the Ml sites, or "pseudo-anions" , to be encapsulated in the same way as the hydrides. 
This also provides a simple physical link between the two LisH structures in the sense that 
they only diflFer in how the chains of hydrides and Ml sites are arranged about the lattice; 
in the Cmc2i structure they are distributed homogeneously while in Abm2 structure they 
are arranged into layers. The chemical difference between the Li^, and the Li^^ atoms is that 
the former bridge two hydrides, whereas the latter bridge two Ml sites. 

In conclusion, evolutionary algorithm searches coupled with first-principles computations 
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predict that a new hydrogen-poor class of compounds, the hthium subhydrides (Li^H with 
m > 1) become stable with respect to decomposition into LiH and Li at P ^ 50 GPa. 
This stabilization occurs at pressures lower than those necessary to form their polyhy- 
dride counterparts. [13j All of the subhydrides are composed of superalkali, LisH, building 
blocks and exhibit ionic bonding on the inside of these clusters and metallic bonding in 
the interstitial regions between them. Analogies were drawn between the geometries and 
electronic structures of the phases predicted in this work and the well-known alkali metal 
suboxides. [TTl [18] The most favorable stoichiometry, LisH, manifested in two iso-enthalpic 
phases which did not exhibit any imaginary phonon modes at 90 GPa. We find the most 
unique aspects of LisH relative to the other Li^H structures we considered to be their ther- 
modynamic and mechanical stability, the mixing of hydridic states with lithium states to 
create subvalence 1-D chain-like electronic behavior, and ordered Li cages around "pseudo- 
anion"-like pockets of electron density resulting in a high symmetry Li sublattice relative to 
the other subhydride structures. Although we have not yet considered finite temperature ef- 
fects nor the kinetic barriers required to break up their would-be precursors (Li or LiH), their 
enthalpic and mechanical stabilities suggest they may potentially be made in high-pressure 
experiments. Furthermore, these results hint that there may be other novel stoichiometrics 
of subhydrides to be explored under pressure and, perhaps, other lithium-rich or suboxide 
materials. 
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